Phase transitions in the frustrated Ising model on the square lattice 
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We consider the thermal phase transition from a paramagnetic to stripe-antiferromagnetic phase 
in the frustrated two-dimensional square-lattice Ising model with competing interactions J 1 < 
(nearest neighbor, ferromagnetic) and Ji > (second neighbor, antiferromagnetic). The striped 
phase breaks a Z4 symmetry and is stabilized at low temperatures for g = J2/I Ji | > 1/2. Despite the 
simplicity of the model, it has proved difficult to precisely determine the order and the universality 
class of the phase transitions. This was done convincingly only recently by Jin et al. [PRL 108, 
045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by 
employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix 
calculations. The J1-J2 model has a line of very weak first-order phase transitions in the whole 
region 1/2 < g < g* , where g* = 0.67 ± 0.01. Thereafter, the transitions from g = g* to g — > do 
are continuous and can be fully mapped, using universality arguments, to the critical line of the 
well known Ashkin- Teller model from its 4-state Potts point to the decoupled Ising limit. We also 
comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin- 
Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity 
of the multicritical point g* in the J1-J2 model. The continuous transitions near g* can therefore 
be mistaken to be first-order transitions, and this realization was the key to understanding the 
paramagnetic-striped transition for the full range of g > 1/2. Most of our results are based on 
Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful 
methodological bench-marks for weakly first-order behaviors and Ashkin- Teller criticality. 

PACS numbers: 64.60. De, 05.70.Ln, 64.60. F-, 75.10.Hk 



I. INTRODUCTION 

The Ising model with nearest-neighbor interactions 
on the two-dimensional (2D) square lattice presents a 
rare instance where the partition function can be com- 
puted exactly at any temperature TX This allows for 
the calculation of the critical exponents characterizing 
the continuous phase transition between the magneti- 
cally ordered ferromagnet and the disordered paramag- 
netic state. Adding competing (frustrated) interactions 
provides a route for the appearance of new phases and, 
in some cases, new types of phase transitions outside the 
Ising universality class. A next-nearest-neighbor anti- 
ferromagnetic interaction represents the simplest way to 
incorporate frustration in the standard Ising model. This 
model, the J1-J2 Ising model, is defined by the Hamilto- 
nian 
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where first and second (diagonal) neighbors on the square 
lattice are denoted by (ij) and ((ij)), respectively, and 
<j 1 = ±1. When the ratio g = J2/IJ1I < 1/2, there is 
an Ising transition versus T to a ferromagnetic stated— 
The competing interactions in the model stabilize a new 
striped phase (see Fig. [I} when g > 1/2. Since these 
stripes can be oriented in cither the x or the y lattice 
direction, the ordering breaks a four-fold (Z4) symmetry 
on the square lattice. Increasing the temperature from 
T = at a fixed g > 1/2, a transition to a disordered 
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FIG. 1. (Color online) The four symmetry related striped 
ground states of the J1-J2 model when g > 1/2. The striped 
phase breaks a Z4 symmetry. Solid and open circles represent 
the spin states Oi = ±1. 

state occurs with no other intermediate broken symme- 
try phase in between. In this paper we study the phase 
transition into the striped state. 

Unlike the Ising transition to a (Z 2 ordered) ferromag- 
netic state, the nature of the phase transition between 
a Z4 ordered state and a disordered state in 2D cannot 
be determined simply from the symmetry of the order 
parameter. This is an example of weak universality, a 
concept first introduced by Suzuki^ where the dimen- 
sionality of the system and the symmetry properties of 
the order parameter are not enough to fix the universal- 
ity and hence, the critical exponents of the phase tran- 
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sition. The exponents may vary with some tuning pa- 
rameter in the system even though the symmetry of the 
order parameter does not change. Only certain ratios of 
the critical exponents remain fixed and these define 7 - the 
weaker form of universality (or, equivalently, the expo- 
nent r\ describing the correlation function at the critical 
point is fixed, while other exponents vary). Some exotic 
2D models, where the critical exponents can be analyt- 
ically calculated as a function of a coupling parameter, 
indeed show this behavior, e.g., the eight-vertex model 8 - 
and the Ashkin-Teller (AT) model^ii 

Even though the frustrated J\-J% model defined by 
Eq. ((T|) represents perhaps the simplest generalized 2D 
Ising model, its stripe transition remained highly contro- 
versial until recently, despite several past studies J 2r&i2r2& 
Early numerical and analytic approaches supported the 
idea that the transition is always continuous for g > 1/2, 
but with critical exponents that vary with g, thus pro- 
viding an example of weak universality. However, some 
variational studie d 13 ' 14 and recent Monte Carlo (MC) 
studie d 15 ' 16 have found a line of first-order transitions 
for 1/2 <g < 1. 

One recent MC study by Kaltz et al.— used the exis- 
tence of a double-peak structure in energy histograms to 
conclude that the transition is first-order up to g = g*, 
with g* ps 0.9. For higher g, in the same work a contin- 
uum field theory was derived perturbatively in 1/g, re- 
sulting in an AT-like model. For intermediate values of g, 
where the (perturbativc) field theory cannot be expected 
to be reliable, and the MC results were ambiguous, it was 
not possible to definitely conclude that the AT scenario 
holds all the way down to g* . In particular, deviations 
from T7 = 1/4 (the fixed value of this exponent in the AT 
model) were seen for g in the range 1 — 5. 

In another recent study^ it was shown by three of 
the present authors that the stripe transition is first 
order in a much smaller range of couplings than pre- 
viously believed; for 1/2 < g < g* , with g* ps 0.67. 
For g > g* it is continuous and in the AT universal- 
ity class. The exponents change continuously with g as 
in the AT models, with g* corresponding to the univer- 
sality of the 4-state Potts mode l 18 ' 19 (which is equivalent 
to the AT model at one end-point of a critical line) and 
g — > oo to standard Ising universality. While AT criti- 
cality had been suspected at the stripe transition earlier, 
it had not been possible to demonstrate this convinc- 
ingly for a wide range of couplings before. The key to 
solving this problem was the realization that the Potts 
model harbors pseudo- first- order behavior and (previ- 
ously known^) logarithmic corrections, and that these 
match very well the properties of the J1-J2 model at 
g ps 0.67. Thus, the full critical curve bridging the Ising 
and 4-state Potts point of the symmetric version of the 
AT model (which we will define in detail further below) 
can be completely realized in the J1-J2 model. 

The pseudo-first-ordcr behavior found in Ref. [l?] im- 
plies that indicators (necessary but not sufficient condi- 
tions) of first-order transitions, e.g., multiple peaks in 



energy and order-parameter distributions, lead to over- 
estimation of the region of discontinuous transitions in 
this model. Mere observation of multi-peak structures 
is not sufficient for concluding that a transition is first- 
order, but careful finite-size scaling studies are required 
to extrapolate, e.g., the latent heat to infinite size. The 
4-state Potts model and neighboring transitions in the 
AT model exhibit clear pseudo-first-order behavior— for 
finite sizes, though these transitions are known to be con- 
tinuous^ It is then necessary to look at certain universal 
properties in the J1-J2 model to determine whether the 
transition is continuous and belongs to the AT univer- 
sality class. This pseudo-critical behavior in the J\-J% 
model was also verified recently in Ref. [2(| where the 
double-peak structure in the energy histogram was shown 
to disappear at large system sizes (L ~ 2000 for a peri- 
odic L x L system) for g = 0.80 (while in Ref. [l7| the 
order-parameter histograms were analyzed). 

In this article we present further evidence to support 
this pictur o 17 ' 20 of the transitions from the striped phase 
in the J1-J2 model. In addition to MC simulations, we 
also consider cluster mean-field (CMF) and numerical 
transefer-matrix (TM) calculations. While in the end 
MC calculations appear to be the only way to reliably 
study the stripe transition close to the most interesting 
point g = g* , due to the subtleties discussed above, it 
is still useful to bench-mark these other commonly used 
methods. 

The rest of the paper is organized in the following 
way: In Sec.|TI]we briefly summarize the known scenarios 
for continuous phase transitions from a Z4 ordered to a 
disordered phase in 2D. We then investigate the phase 
transitions of the J1-J2 model in detail using the cluster 
mean- field theory approach (Sec. IIIip . MC simulations 
fSec.HVp and the TM approach (Sec.|V]). We also present 
some further results for the AT model in Sec. IIV1 includ- 
ing its pseudo- first-order behaviour near the 4-state Potts 
point. We further establish the equivalence between the 
continuous phase transitions in the AT and J1-J2 models, 
including quantitative results for how the parameters of 
the two models correspond to each other in terms of the 
varying critical indices. We give a brief summary of the 
results in Sec. IVI1 



II. EXPECTATIONS FROM UNIVERSALITY 

In two dimensions, the critical exponents can have var- 
ious possible values when the ordered phase breaks a Z± 
symmetry. In the J1-J2 model for g > 1/2, only the 
g — >• 00 limit and 5 = 1/2 transitions are exactly known. 
At g —¥ 00, the system consists of two decoupled Ising 
systems and there is a continuous thermal phase transi- 
tion in the Ising universality class. At g = 1/2, it is clear 
that there is a first-order transition at T = between 
a ferromagnetic and a stripe-antifcrromagnetic state of 
the type depicted in Fig. [T] The first-order transition 
point is unusual in that there is a co-existence of a large 
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number of states^ made up entirely of horizontal (or, 
vertical) stripes where the orientation (a = +1 or —1) of 
each stripe can be chosen independently The nature of 
the g > 1/2, T > transitions is not a priori clear. We 
briefly discuss two microscopic scenarios which cover the 
known theoretical possibilities for continuous transitions 
out of a Zi symmetry-broken 2D state. 

Let us first consider the 2D XY model in a four-fold 
anisotropic field of strength hi, 

H = - cos ( e * h * /Z cos ( 4 ^)' ( 2 ) 

where the sites i reside on a square lattice and 9i de- 
fines a 2D fixed-length vector in the XY plane. At 
hi = 0, there is a Kosterlitz-Thouless (KT) transition 
versus temperatur o 21 i 22 while \h4,\ — > oo gives the stan- 
dard Ising universality. A non-zero /14 leads to a four-fold 
broken symmetry phase at low T. The critical exponents 
change as a function of hi, e.g., the thermal exponent v 
equals 1 in the Ising limit and v — > 00 at the KT transi- 
tion. However, the important observation for our purpose 
here is that the specific heat exponent a/v pertaining to 
finite-size scaling equals in the Ising limit (the specific 
heat diverges logarithmically with system size here) and 
develops a cusp at finite /14 indicating a negative a. This 
cannot possibly explain the behavior of the specific heat 
in the J\-Ji modeUi where the divergence with system 
size seems quite strong for all g > 1/2, indicating a/v > 
if the transition is assumed to be continuous. 

Next, we consider the AT model on the square lat- 
tice2»ii which can be written as 

H = - y^(cr,-<7j + TiTj + l\n,a,T,r, ). (3) 

where two Ising variables a% , Tj reside on each site i of the 
square lattice and are coupled to each other through K. 
There is a symmetry of the model, corresponding to the 
permutations of the variables a, r, and o~t. These map 
the Hamiltonian © onto itself, and, thus, only values of 
K in the range [—1,1] have to be considered. 

The ferromagnetic phase of the AT model breaks a Z4 
symmetry and is defined by (err) 7^ and (a) = ±(t). 
The transition from this Z± ordered state to the fully dis- 
ordered state ((err) = and (a) = (r) = 0) has continu- 
ously changing exponents which are exactly known as a 
function of K 10 ^ 11 using the following relations based on 
the powerful Coulomb-gas (CG) formulation for study- 
ing this class of 2D phase transitions (see Ref. [h] for an 
excellent review of this approach): 

y t =2-2/g R , j, h = 15/8, y p = 2 - l/(2g R ), (4) 

with 

g R = ® arcsin Q coth(2/T c )^ (5) 

being the CG coupling. The critical temperature T c is 
exactly given by the fact that there is a self-dual line; 

sinh(2/T c ) = cxp (-2K/T C ). (6) 



The three exponents are the thermal exponent y t = 1/v, 
the magnetic exponent y^ = 2 — r)/2, which is fixed in 
the region, and the exponent y p corresponding to a po- 
larization field acting on one of the sets of Ising vari- 
ables, PJ2i T i (which breaks the Z2 symmetry of the 
Hamiltonian).— The corresponding scaling dimensions 
are X t = 2 - y t , X h = 2 - y h , X p = 2 - y p . 

In the AT model, K = corresponds to the decou- 
pled Ising limit and K = 1 corresponds to 4-state Potts 
model universality. When K is extended from to neg- 
ative values, the thermal exponent v increases and the 
specific heat develops a cusp. In the AT model, the crit- 
ical line is defined up to the point K = — 1 where v = 2. 
The transitions from K = to K = —1 can be viewed 
as a subset of the critical points in the anisotropic XY 
model discussed above. Since the exponent a/v > 
when K € [0, 1], this suggests that if there are contin- 
uous phase transitions in the J1-J2 model, then these 
belong in the same universality class as critical points in 
the AT model in the range K £ [0, 1]. Moreover, the spe- 
cific heat exponent a/v seems to decrease smoothly as g 
is increased^ and its behavior with g indicates that the 
frustrated Ising model may have all the critical points 
of the AT model from K = till K = 1 (4-state Potts 
model), which will then be the multicritical point in the 
frustrated Ising model and is the end point of the line 
of continuous transitions between Z4 ordered and disor- 
dered phases in the AT model. 

As we will see later in Sec. IIV1 the continuous tran- 
sitions in the J1-J2 model for g £ [g*,oo) indeed can 
be mapped to the transitions in the AT model when 
K e [1,0). We will asume that the stripe transition 
for g < g* is first order instead of some unlikely and al- 
ternative exotic behavior outside the known scenarios for 
Z^ symmetry breaking. Results near g = 1/2 discussed 
in Section [TV] suggest that the first-order transitions are 
very weak, and large system sizes may be needed to ob- 
serve unambiguously the expected first-order scaling be- 
haviors. In principle one cannot exclude, based on the 
numerics alone, that there is some yet unknown type of 
continuous transition for 1/2 < g < g*. 

III. CLUSTER MEAN-FIELD THEORY 

Here we study the J1-J2 model using a CMF approach 
based on 2 x 2 and 4x4 clusters, as illustrated for the 
latter case in Fig. [5] In this section we will not just 
consider the striped phase and transition, but also the 
standard ferromagnetic phase obtaining for < g < 1/2. 

A. Variational approach with a reference system 

One way to formulate a mean-field theory is to con- 
struct an approximate expression for the partition func- 
tion with the aid of some solvable models Let E® be the 
energy for such a reference system in spin configuration 
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FIG. 2. (Color online) Illustration of a variational mean-field 
theory based on 4 x 4 clusters. The infinite lattice is divided 
into clusters (sites connected by the thicker, red lines) , and a 
single isolated cluster with added magnetic fields (indicated 
here by different circles) is solved exactly (with no interactions 
between it and neighboring clusters), giving (E°)o and F°. 
The energy {E c )o of the original system without the fields 
is evaluated using the mean-field decomposition ((Tiff j}o — > 
{ci)o((Jj}o for the bonds connecting clusters and using the 
imposed periodicity to translate both sites i and j into the 
same cluster. The fields are adjusted to minimize the upper 
bound F* — F® + (E® — E c )o on the cluster free energy. 



a. It is assumed that its partition function, 



^o = E 



(7) 



can be calculated exactly in some way (numerically or 
analytically) . We can write the partition function of the 
original system of interest, with energy function E ai as 



z = z J2 



-E°/T 



Zo(e -( S -^)/T )oj 



(8) 



where ()o denotes the expectation value with respect to 
the Boltzmann distribution of the model E®. Under our 
assumption, this expectation value can also be evaluated 
exactly A well known result is that22 



(9) 



from which we obtain an upper bound F* to the free 
energy F = -Tln(Z); 



F < F* = F + (E — E°} . 



(10) 



This variational principle for the free energy is very useful 
if the reference model E® has some parameters that can 
be varied. One can then minimize the free-energy bound 
F* with respect to those parameters, to obtain the best 
(in the sense of minimum free energy) variational solution 
to the system E a . 

In the variational CMF approach the reference system 
is a infinite system is divided into clusters, with no inter- 
actions between the clusters. The energy of the infinite 
reference system can be written as a sum over the iden- 
tical clusters c; 



E [) 



(ii) 



A small isolated cluster can be solved by exact summa- 
tion over all its spin configurations. The aim is to mini- 
mize the free energy of the J1-J2 model with respect to 
variational parameters of the reference system. In prin- 
ciple, the reference model can contain arbitrary field and 
spin dependent terms within the cluster. However, in 
practice, the cluster Hamiltonian function minimizing the 
variational free energy has exactly the same couplings as 
the original J\-Ji model, and only fields —hi<Ji acting on 
the edge spins are added. Since the J1-J2 model includes 
only up to second-neighbor couplings, the edge here has 
the standard meaning of only the outermost layer of sites, 
but for longer-range interactions the "edge" extends fur- 
ther into the cluster. 

Let us first consider independent fields hi coupling to 
all the spins <7j , i = 1 , . . . , n within a cluster of n sites 
(disregarding symmetries that will eventually imply that 
some of the fields should be equal); 



E° = JijWj ~ E h * a * 

(ij) «'=1 



(12) 



where (i, j) refers to site pairs (counted only once) within 
the cluster. In the model we will consider explicitly, 



-J\ or J2, but in principle Jij could include even 



longer-range interactions. The cluster energy defines the 
reference Boltzmann distribution with relative probabil- 
ities Wo(ct) = exp[— E°(a)/T] for the 2™ spin configura- 
tions a = <7i, . . . , a n . Using this probability distribution, 
we can evaluate the partition function Z® = Wq(<j), 
(Eq)q, and the expectation value (E)q of the original en- 
ergy for each cluster. We have 



(E c )o = E Jij(<7i°j)o + \ Jij{<Ti)o(<7j)o, 



(13) 



(ijy 



where in the second sum refers to interactions be- 

tween a site i in the cluster c and a site in a different 
cluster. Since all clusters are equivalent, this site can 
be translated into an equivalent site j of the cluster c. 
The factor 1/2 accounts for the fact that each interac- 
tion bond is shared by two different clusters. 

In the free-energy bound {TDJ F* = -Tln(Z°) + (E - 
E°)o we need only the difference between the two ener- 
gies, for which the intra-clustcr interactions cancel; 



(E - E°) Q = ^Y1 J v( a M°j)o + zZ M^)o- 



(14) 
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To minimize F* we need its derivatives with respect to 
the fields; 



Q F * 



Jij(<7i)o(((7jcrk)o - (o-j)o(o-fc)o), (15) 



(ijy 
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FIG. 3. Clusters and fields used in mean-field calculations for 
the J1-J2 Ising model. For a ferromagnetic state (all black cir- 
cles), all fields are positive, of strength hi for the 2x2 cluster 
and h\,h,2 for the 4x4 cluster (with ho — 0, because the 
sites marked are not on the cluster edge). For the striped 
state, the black and white circles indicate positive and nega- 
tive fields of magnitude h s and sign (— l) Xi , where Xi is the 
^-coordinate of the sites i to which it couples. For the 4x4 
cluster, the broken rotational symmetry of the striped state 
implies one more variable field than for the ferromagnet, for 
a total of three adjustable parameters hi, hi, /13. 



which can be written in the form 

OF* / \ 

g^=J2{ hi + J2 J ^( (T ^)a ik = 0, (16) 

where the notation in the second sum sum indi- 
cates summation for given i over only those spins j 
corresponding to inter-cluster (edge) interactions and 
<Hk = (ciCfc)o — ( cr i)o( cr fc)o- These equations are satis- 
fied if 

hi = J ij{vj}o, (17) 

Uh 

which amounts to self-consistency conditions for all the 
fields. For sites i that have no non-zero inter-cluster in- 
teraction Jij, we have hi = 0, i.e., we need to consider 
only fields on the edge spins. This self-consistent solution 
has the lowest free energy also if other interactions are 
allowed within the reference system E°. The variational 
approach is therefore equivalent to the self-consistent ap- 
proach. The advantage of starting from the variational 
ansatz is that the free energy (its upper bound F*) is 
also obtained without any further assumptions. 

In practice one does not have to treat all the fields hi 
as independent parameters, because the optimal fields for 
an ordered state will obey symmetries corresponding to 
those of the order parameter. For the J\-Ji Ising model 
considered here, we have ferromagnetic and striped or- 
der for J2/J1 < 1/2 and > 1/2, respectively. The field 
arrangements appropriate for these order parameters on 
2x2 and 4x4 clusters are illustrated in Fig. [3J We will not 
consider the 3x3 cluster, because it is not appropriate for 
the striped state (due to its incompatibility with the pe- 
riodicity 2 in one of the directions, although in principle 
one could also take this into account by modified bound- 
ary conditions). Using the appropriate symmetries, the 
largest number of parameters is here three, for the stripe 
order on the 4x4 lattice. With such a small number 
of parameters, we can easily solve the self-consistency 
equations numerically. 
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FIG. 4. (Color online) Phase diagram of the J1-J2 Ising model 
in the coupling-temperature plane based on mean-field calcu- 
lations with clusters of size 2x2 (thin curves, at higher T) 
and 4x4 (thick curves, at lower T). Dashed and solid curves 
indicate continuous and first-order phase transitions, respec- 
tively, with the circles indicating the multi-critical points at 
which the order of the transition changes. 



B. Phase diagram 

By finding the optimal solutions for both ferromag- 
netic and striped field patterns, and comparing their 
free-energy minimums for a range of temperatures and 
coupling rations <?, the phase diagram of the system can 
be extracted. To precisely determine a phase boundary 
as a function of temperature at fixed g, one can carry 
out a bracketing procedure to locate the point at which 
the optimal solution changes between paramagnetic and 
ordered, or between ferromagnetic and striped ordered. 
Fig. [4] shows the phase diagram obtained this way, based 
on both 2x2 and 4x4 clusters. There are continuous 
as well as first-order transitions. First-order transitions 
obtain close to g = 1/2, which is the point at which we 
already concluded that there should be such a transition 
as a function of g when T — > 0. The point at which the 
transition becomes continuous is stable with respect to 
the cluster size, g* « 0.66, and is in remarkably good 
agreement with the value g* ~ 0.67 obtained in the pre- 
vious MC workii identifying the Potts point. The mean- 
field calculation can of course not give any information 
on the true critical exponents. 

The paramagnetic-ferromagnetic transition is seen to 
always be continuous within the 2x2 cluster calculations, 
but it also changes to first-order in a narrow window of 
g values when the 4x4 cluster is used. There are no 
clear indications from previous MC simulations of the 
transition being first-order in this regime, but it may be 
worth examining this issue more carefully as well with 
improved MC simulations. In this paper we focus on the 
stripe phase for g > 1/2, however. 

In principle one can try to extrapolate the critical tem- 
perature to infinite cluster size, but this is not possible 
based on just the 2x2 and 4x4 results. In principle 
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FIG. 5. (Color online) Self-consistent variational field param- 
eters (left) and the corresponding induced stripe magnetiza- 
tions (right) versus temperature obtained using a 4 x 4 cluster 
at J2/J1 = 0.7. The inset of the right panel shows the ratio 
of the shell magnetizations for shells s > to the one with 
s = (center of the cluster). The shell index s follows the 
convention illustrated in Fig. [3] 
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FIG. 6. (Color online) Free energy versus the external field 
parameter h in the 2x2 cluster mean-field theory at and close 
to first-order transitions. Positive and negative h correspond 
to ferromagnetic and striped field patterns, respectively, (a) 
shows results for different coupling ratios at and close to 
the ferromagnetic-striped transition at fixed T/Ji — 1, and 
(b) shows the behavior for temperatures at and close to the 
striped-paramagnetic transition at g = 0.55. 



one can add standard single-site MF calculations and the 
3x3 cluster to enable some estimates. Here, however, 
our interest is in the order of the transitions and we will 
not attempt any CMF-bascd extrapolations of T c (which 
anyway come out very precisely in the MC and TM cal- 
culations, as discussed in the later sections). 



C. Order parameter and free energy 

An example of self-consistent field parameters and in- 
duced shell magnetizations is shown in Fig. [5] These 
results are for the 4x4 cluster at g = 0.7, where 
there is a continuous paramagnetic-striped transition at 
T/Ji k, 1.83. All the fields vanish continuously at this 
point. Note that the order parameter is not uniform, as 
it should be for an infinite system, but shows significant 
variations between the shells. For most of the tempera- 
ture range, the order is the weakest at the four central 
spins, where there is no field, but close to T c one of the 
edge magnetizations becomes equal to it. The ratio of 
the magnetizations m s to the central magnetization mo 
is significantly different from 1 close to T c , but should re- 
main finite because all m s should vanish as T — > T c with 
the mean-field power-law behavior m s ~ (T c — T) 1 / 2 . 
In general, one would expect that the four central spins 
should most closely represent the behavior of the infinite 
system. If we could go to much larger cluster sizes, we 
would expect the order parameter to become uniform in 
the interior of the cluster, with non-uniformity emerging 
gradually as the edges arc approached. 

To discuss the first-order transitions, it is useful to ex- 
amine the free energy of the 2x2 cluster, where there is 
just one variational parameter. The free energy for both 
the ferromagnetic and striped fields can can be showed 
in the same graph by defining a new parameter h, such 
that for ft > this is the ferromagnetic field hi = h, 
whereas for ft, < it is the strength of the stripe field; 
hi = \h\. Two examples of the dependence of the cluster 
free energy F* on this parameter as a first-order tran- 
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FIG. 7. (Color online) Dependence on the coupling ratio of 
the discontinuity of the ferromagnetic magnetization (left) 
and stripe magnetization (right) at the first-order transitions 
obtained with the 2x2 and 4x4 clusters (I — 2, 4). Note in the 
left panel and in Fig. [4] that for g in the range 0.50 ~ 0.504, 
as T is lowered there is afirst a paramagnetic-ferromagnetic 
transition, followed by a transition into the low-temperature 
stripe state. 



sition is crossed are shown in Fig. [6] In the left panel, 
two minimums for h ^ can bee seen, corresponding 
to ferromagnetic and stripe orders, whereas in the right 
panel one of the minimums is at ft = 0, corresponding 
to the paramagnetic phase, and the other minimum is 
for a striped state. In either case, the two minimums 
are degenerate at the transition between the two phases. 
Changing g at fixed T (as in the left panel) or varying T 
at fixed g (right panel) , the degeneracy is broken and one 
of the states becomes the stable one. The other, higher 
minimum then corresponds to a meta-stable state. 

The discontinuities associated with the first-order tran- 
sitions vanish continuously at the special multi-critical 
points indicated with circles in the phase diagram in 
Fig. |U Fig. [7]shows the behavior of all the discontinuities 
for the 2x2 and 4x4 clusters. 
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FIG. 8. (Color online) (a) Universal crossing value U* of the 
Binder cumulant vs K for the AT model and vs l/(l + <? — <?*), 
g* — 0.67 for the J1-J2 model. The AT data points have 
been fitted to a single polynomial function, while the J1-J2 
data were fitted with several polynomials in segments. The 
corresponding K values with the same U* as g = 1.0 and 
g = 2.0 are marked in the graph. The map between the two 
models at these points are: (g — 1.0) l/(l + g — g*) = 0.752 to 
K = 0.35 with U* = 0.8805; (g = 2.0) 1/(1 +g- g*) = 0.429 
to K = 0.081 with U* = 0.9086. (b) Map between the AT 
model (K) and the J1-J2 model (g~ 1 ) using the universality 
of U* and the procedure illustrated in (a). 



A. Calculated Observables 

Before proceeding further, we define the observables 
that we measure in our MC simulations for the J1-J2 
and the AT models respectively. Let us first define the 
order parameters that characterize the broken Z4 phase 
in both the models. The striped phase of the 2D frus- 
trated Ising model is characterized by a two-component 
order parameter (m x ,m y ) with 



i=l 
1 N 



(18) 
(19) 



i=i 



where (xi,yi) are the coordinates of site i on a L x L 
periodic square lattice and TV = L 2 . We define to 2 = 



to 2 + to 2 and the stripe susceptibility as 

x = y «™ 2 >-<H> 2 )- 

We also measure the specific heat, 

C v = ^((E 2 ) - (E) 2 ), 

where E is the energy per site. For the AT model, the 
order parameter can again be expressed as a 2D vector 
(m CT , m T ) where 



(20) 



(21) 



i=l 
1 N 



(22) 
(23) 



IV. MONTE CARLO SIMULATIONS 



We have simulated the J1-J2 Ising model using a stan- 
dard single-spin Metropolis algorithm^ Due to the pres- 
ence of frustration, cluster MC methods cannot be used 
for this model unless J\(J2) = 0. We found that single- 
spin Metropolis algorithm is sufficient to study the ther- 
mal phase transitions accurately if g is not very close to 
1/2 (we have gone up to g = 0.52 using single spin-flip 
MC moves). The transitions closer to g = 1/2 can be 
simulated using a combination of parallel tempering and 
certain non-local spin flips^ which have a high accep- 
tance probability very close to g = 1/2. We have also 
simulated the AT model on the square lattice, and for 
that we again use a local Metropolis algorithm except at 
K = 1, where we use a cluster algorithm^ Temperature 
is measured in units of J\ for the frustrated Ising model 
and in units of K for the AT model. MC simulations 
combined with finite-size scaling and universality argu- 
ments provides the most unbiased method to understand 
the nature of the transitions in the 3\-Ji model. 



and m 2 = m 2 , + m 2 .. The order parameter susceptibility 
X and specific heat C v are then defined exactly in the 
same way as described above for the J1-J2 model. 

We also compute the Binder cumulant of the order 
parameter in both models. It is defined as 



(7 = 2 1 



1 (to 4 ) 

2 (to 2 ) 2 



(24) 



where the constants are chosen to give a step-function 
(U — > and U — >• 1 in the disordered and ordered phase, 
respectively) for a 2D vector order parameter in the ther- 
modynamic limits 

Lastly, we collect the histograms of the squared order 
parameter m 2 and the energy E near the transition for 
both the models. These are helpful for analyzing the 
pseudo-first-order behavior in detail. 



B. Map between AT and J1-J2 critical points 



The Binder cumulant, Eq. (|24|) . turns out to be es- 
pecially useful in establishing the universality class of 



8 



10' 



10" 



C max 










' — g=2.0 




o-o K=0.35 




«g=1.0 - 




10 20 



50 100 200 



0.01 



10 20 



50 100 200 



FIG. 9. (Color online) Examples of divergent peak-values of 
the specific heat and stripe susceptibility vs L for the J1-J2 
and AT models. The factor L 7//4 corresponding the asymp- 
totic AT scaling of the susceptibility have been divided out 
in the right panel. Two point pairs are chosen for which the 
mapping in Fig. [5] show that the critical exponents for the two 
models should coincide. 



the continuous transitions in the J1-J2 model. It is 
well knowr>2i that for continuous phase transitions, the 
Binder cumulant for different system sizes cross at the 
critical point (when the system size is large enough). The 
value of the crossing U* is universal as well and charac- 
terizes the universality class of the phase transition. U* 
may in some cases depend on details of the model beyond 
the universality class, like the boundary conditions and 
shape of the lattice and the anisotropy of the interac- 
tions^ However, in our case, both the J1-J2 model and 
the AT model live on periodic square lattices and the in- 
teractions respect the full symmetry of the lattice, so a 
comparison of U* between the two models seems to be 
justified. We have established this directly from our MC 
data by using the equality of U* to map phase transi- 
tions in one model to the other and then directly looking 
at critical exponents to check if the universality class is 
indeed the same. 

We estimate U* from our MC simulations by extract- 
ing the crossing point of U between data for pairs (L, 2L) 
and then extrapolating to L — > In Fig. [8] we show 

[/* as a function of the coupling K and g for the AT 
model and the J1-J2 model, respectively. This immedi- 
ately establishes a numerical map between the continu- 
ous transitions of both the models. From Fig. [8j we see 
that g 0.67 corresponds to the 4-state Potts model 
universality in the AT model (K =1). This was already 
reported in Ref. [I?] and other consistency checks were 
used there to show that the multicritical point is located 
at g* = 0.67±0.01. Since then, the location of g* has also 
been computed in Ref. H3, and the result agrees perfectly 
with the earlier result. 

As a further illustration of the correctness of the proce- 
dure, we use the data presented in Fig.[5]to note that the 
phase transition at g = 1 should map to K sa 0.35 and 
g = 2 to K « 0.081. This is indeed consistent with the di- 
vergent critical forms of physical quantities. In Fig. [5] we 
plot the peak value of the specific heat C max (£) and the 
order parameter susceptibility Xmax(^) versus L for the 



two models at the above points. By standard finite-size 
scaling arguments, C max - L a l v and Xmax(i) ~ £ 7/t '■ 
For the system sizes studied here (L < 256), the expo- 
nent a/v for the two models, estimated from the slope of 
C max (L) on a log- log scale (Fig. [9l left panel), converges 
to the same value in a very similar way for the J1-J2 
and Potts model at the corresponding g and K values. 
This is also the case for the exponent j/v (Fig. HO right 
panel), which converges to the value 7/4, as is expected 
for AT universality. The latter behavior is of course less 
useful, since the exponent 7 = 7/4 is expected on the 
whole critical curve. In both cases some deviations from 
pure power-laws can be observed, and we will discuss this 
further below. The behavior seen for the specific heat is 
nevertheless quite telling and suggestive of the same crit- 
ical exponent in the two models at the mapped points. 
Thus, we have rather convincingly established the map 
(Fig. [SJ between the parameters of the AT model and the 
frustrated 3\-Ji Ising model in two different but mutu- 
ally consistent ways. Note that this kind of map does not 
imply microscopic equivalence of the two models, which 
holds only in the weakly coupled Ising limit (l/p — ► 0), 
but demonstrate common low-energy descriptions of the 
systems for the mapped parameter values. 



C. 4-state Potts scaling at g* 

To further strengthen the case for<?*=0.67±0.01 be- 
ing in the 4-state Potts universality class, we next con- 
sider scaling in the form of data collapse of the specific 
heat and the susceptibility at the coupling g = 0.68 of the 
J1-J2 model (for which we have MC data; the estimate 
for g* is based on interpolation, as shown in Fig. [8j and 
for the 4-state Potts model on the square lattice (K = 1 
in the AT model) . The critical exponents of the 4-state 
Potts model areJ£ v = 2/3, a/v = 1,7/1/ = 7/4. There 
are however important multiplicative logarithmic scaling 
corrections at this critical point that strongly affect finite- 
size scaling. According to Ref. Qjl the divergences in the 
thermodynamic limit are of the forms: 



X 



e~ir 2/3 (-iog|ii) 1/2 



(log0 3/2: 

£7/4 

(log£) 1/8: 



(25) 

(26) 
(27) 



where t = (T — T c )/T c is the reduced temperature 
and £ is the correlation length. Then, using finite-size 
scaling arguments, one expects C v L~ 1 (log(L/L ))~ 3 / 2 
and xL~ 7 / 4 (\og(L/Lo)) 1 / 8 to be functions of the ar- 
gument t(— log |t|) -3 / 4 L 3 / 2 . The two quantities with 
these expected asymptotic L-dependences divided out 
are graphed versus £(— log |i|)~ 3 / 4 L 3 / 2 for different sys- 
tem sizes L in Fig. 1101 with the non- universal scale factor 
L (which can be different for different quantities) treated 
as a fitting parameter to optimize the data collapse. 
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FIG. 10. (Color online) Data collapse with the anticipated 
leading logarithmic correction to the specific heat C and the 
susceptibility % f° r the 4-state Potts model and the J1-J2 
model at g = 0.68. The system sizes included in (a) are 
L = 60, 64, 80, 120, 128, 160, and in (c) L = 48, 64, 96, 
128, 192, 256, while for (b) and (d) all the sizes are listed in 
the panels. The black curves are common curves fitted to all 
data points shown in each figure. The corresponding reduced 
chi-square per degree of freedom for the fitted curves in (a), 
(b), (c) and (d) are X 2 = 1-3, 1.6, 1.2, 1.8. 



.L=128 (a) 
\P(m) mL =256 
L=384 
L=512 





-1.60 -1.55 -1.50 -1.45 -1.40 



14 
12 


: p(E) 




10 






8 
6 
4 






2 







-1.10 


-1.05 -1.00 -0.95 
E 



FIG. 11. (Color online) Histograms of the squared order pa- 
rameter m 2 and the energy E for (a), (b) the 4-state Potts 
model (the K = 1 AT model) and (c),(d) for the J1-J2 model 
at g — 0.67. Here T is very close to T c , chosen such that the 
two peaks in the energy histograms are of the same height; 
for the Potts model T/K = 3.64231 for L = 128 and 3.64460 
for L = 256, while for the J1-J2 model T/Jl = 1.2014 for 
L = 128 and 1.2004 for L = 256. 



Fig. HUla) shows the resulting data collapse of C v for 
the 4-state Potts model with a set of moderate to large 
system sizes L = 60 — 160 included in the fitting proce- 
dure (with smaller sizes excluded because they are visibly 
affected by subleading corrections to scaling). The data 
collapse to a common fitted polynomial is statistically 
sound with the parameter L = 0.20 ± 0.01. Fig. flDTb) 
shows the same kind of analysis for the J1-J2 model at 
g = 0.68. The system sizes included here are in the range 
L = 80 - 128 and L = 0.144 ± 0.006. The data collapse 
of C v in both the cases, using the same expected expo- 
nents and multiplicative logarithmic corrections, is an- 
other strong indication that g = 0.68 in the J1-J2 model 
is in close neighborhood of the 4-state Potts end-point of 
the critical Potts-Ising line of the AT model. 

The logarithmic scaling correction for the susceptibil- 
ity, Eq. (|2"T|) , does not yield a good data collapse for either 
g = 0.68 or the 4-state Potts model. Therefore, instead 
of using xL~ 7 / 4 {\og L/Lq) 1 / 8 on the y-axis, we treat the 
exponent of the logarithmic function as another variable 
r in addition to Lq. After carrying out a multi- variable 
data collapse, we obtain r close to —1/8, instead of the 
value 1/8 proposed in Ref . [l9|. Figs. [TUtc). (d) show the 
results for the two models with r = —1/8. Given that the 
data collapse is very good in both cases, and only the sign 
of the log-exponent differs from what was expected, the 
natural conclusion is that there is a sign mistake in the 
analytical result for this exponent in in Ref. [l9l Thus, 
we propose that \ ~ £ 7/,4 (l°g£) 1/ ' 8 ! instead of the form 
in Eq. dUD. 



D. Pseudo-first-order behavior in the AT model 

Even though the transitions in the AT model are all 
continuous, there are interesting pseudo-first-order signa- 
tures at finite system sizes at the 4-state Potts point and 
its neighborhood. This was explicitly shown in Rehous- 
ing various observablcs. Here, to further investigate this 
behavior, we show histograms of the distribution of the 
squared order parameter m 2 and energy per site E for the 
4-state Potts model in the top panels of Fig. [TT] For each 
system size, the temperature is very close to T c , chosen to 
ensure that the two peaks in the energy histograms are of 
the same height. There is clearly a double-peak structure 
present even for the large system sizes L = 512. The dis- 
tance between the peaks decreases slowly as the system 
size increases and the dip between the peaks also does 
not increase appreciably (as it should if the double-peak 
structure is evolving into delta-functions corresponding 
to phase coexistence at a first-order transition). This can 
in principle happen for a weak first order transitions as 
well, if system sizes used are L <C £, where £ is the large 
but finite correlation length at the transition. However, 
in this case we know rigorously that the 4-state Potts 
model harbors a continuous transition. For a continu- 
ous transition, there cannot be an order parameter jump 
or a latent heat in the thermodynamic limit. Thus, the 
distance between the double peaks will eventually shrink 
to zero when L — > 00. This type of double peak struc- 
ture was previously also observed in the energy distribu- 
tion of the Baxter- Wu model on the triangular lattice^ 
which is in the same universality class as the 4-state Potts 
model. The corresponding histograms of the J1-J2 model 
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at g = 0.67 (which equals g* within error bars) also show 
a very similar behavior as can be seen from figures in 
the bottom panel of Fig. [TTJ again confirming the equiv- 
alence of the critical behaviors with that of the 4-state 
Potts model. 

In Ref. [I?], it was also shown that the Binder cumulant 
of the order parameter exhibits a non-monotonic behav- 
ior with T, developing a negative peak at K = 1.0 and 
its vicinity (e.g., at K = 0.95) in the AT model. A nega- 
tive Binder peak is often taken as evidence of first-order 
transition, but here the transitions are clearly continu- 
ous. However, the negative peak increases very weakly 
with system size L, with the increase being much slower 
than the expected^ L 2 divergence. Moreover, the dip 
is more pronounced at K = 1 compared to K = 0.95, 
which indicates that there may be a K* below which 
these pscudo-first-ordcr signatures vanish. These pscudo 
first order signatures in turn lead to an over-estimation of 
the region of first-order transitions when the appearance 
of double-peak structures in energy or order-parameter 
histograms are taken to be indicative of discontinuous 
transitions in the J\-Ji model (see Ref. [l?] for more dis- 
cussion on this point). 

Note that a two-peak structure was found in the energy 
histogram at very large system sizes for g = 0.9 in Ref.ll6l. 
which was taken as evidence of the first-order transition 
extending at least up to this value; based on our con- 
clusions we instead see that the pseudo-first-ordcr region 
in the J1-J2 model extends from g* « 0.67 to g < 1. 
Following our previous work, Ref. [20| recently considered 
the energy histograms at g = 0.8 for even larger system 
sizes than before. A double-peak structure appears when 
the system size reaches about L = 1000, but going fur- 
ther it eventually disappears again, around I = 2000. 
This again confirms the pseudo-first-order behavior the 
J1-J2 model close to g*, and is a demonstration of pit- 
falls in distinguishing between continuous and first-order 
transitions. We stress here that the reason we were able 
to avoid this pitfall is that the Potts model is rigorously 
known to harbor a continuous transition, and we found 
that the behavior of the J1-J2 model at g — g* « 0.67 
matches it very well in all respects. Thus, the combina- 
tion of analytical theory and numerics was crucial. 



E. Weak first-order transitions 

Since the J\-3i model at g* is in the 4-state Potts uni- 
versality class and approaches the Ising limit for higher 
<7, the transitions for 1/2 < g < g* have to be first-order 
transitions if there is to be a correspondence with the 
known scenarios for ^-breaking transitions discussed in 
Sec. [TTJ Unless some unknown scenario applies, which we 
find unlikely (but cannot rule out completely), all the 
transitions in the range 1/2 < g < g* are very weak first- 
order transitions. Weakening of a discontinuous transi- 
tion is expected when approaching a multicritical point, 
since the continuous transition has to be approached in 
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FIG. 12. (Color online) Scaling exponents a/u and j/u of 
C m ax(i) and Xmax(i). The exponent a/v is calculated based 
on the local slope of C max (L) between L and L/2; j/v is 
calculated similarly from Xmax (L). The system sizes L < 128 
in (a) and (c); and L < 256 in (b) and (d). 

a continuous manner. However, here the transitions in 
the close neighborhood of the obvious first order point 
g = 1/2 arc also weakly first order. Ref. [ID suggested 
this based on the appearance of double peak structure in 
the energy histograms. However, as we saw, such double- 
peak structure also appears at the Potts point and its 
neighborhood (where they disappear in the thermody- 
namic limit). Here we show the evolution of the pseudo 
critical exponents with system size L on the (likely) first- 
order side close to the g — 1/2 point. 

We analyze the peak value C max (L) of the specific heat 
and Xmax (I) of the stripe susceptibility. By finite-size 
scaling arguments for first-order transitions, these quan- 
tities should diverge as L 2 in 2D.— Examples of the scal- 
ing behavior are shown in Fig. Q21 Two coupling ratios 
g = 0.52 with system size L < 128 and g = 0.55 with 
system size L < 256 are considered. The peak value 
of C max (I) and Xmax (L) are shown as the inset in each 
graph on a log-log scale. Graphed in this way, the peak 
value of C max (I) and Xmax(L) seem to follow a linear 
scaling behavior, especially in the insets of Figs. [TSTb), 
(c), (d). A more systematic analysis involves extract- 
ing the running exponents ^(1) and ^(1) from the local 
slope of C max (I) and Xmax(I) between, e.g., system sizes 
L and L/2. These should approach 2 as L — > 00 for a 
first-order transition. The first-order exponent 2 is not 
obtained in these figures for the system sizes studied. 
The scaling exponents f(L) and ^(L) [Figs. IT2T a). (c)] 
increase as the system size increases for g = 0.52, but 
have not converged at size L — 128. This suggests that 
it may require very large system sizes to observe the ex- 
pected L 2 scaling behavior. The scaling exponent ^(1) 
for g = 0.55 [Fig. IT2T b')] is further away from 2 at the 
same system size L = 128, while it shows the same ten- 
dency to increase as g = 0.52. The scaling exponent -(I) 
for g = 0.55 actually seems to have in fact converged to 
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w 1.82, at a comparably smaller system size L = 80. This 
is again may be indicative of the large correlation length 
involved at a very weak first-order transitions, this being 
reminiscent of extremely weak first-order transitions like 
that in the 5-state Potts model^i on the square lattice. 
Itt is puzzling, however, that there is no visible size de- 
pendence in the local -^{L) between L = 80 and 256 for 
g = 0.55. It would be interesting to go to larger lattices. 

Note that for both couplings — (L) is already much 
larger than 1, which is the maximum divergence expected 
if the critical point is continuous and in the AT universal- 
ity class. It is an interesting open question to understand 
the mechanism which makes the transitions near g — 1/2 
so weakly first-order; it likely is related to the extensive 
degeneracy of the system at T = 0. 



V. TRANSFER MATRIX CALCULATIONS 

We now address the stripe transition in the 3\-Ji 
model using numerical TM calculations. Consider the 
lattice wrapped on a cylinder of infinite length with cir- 
cumference L, on which the TM is constructed. It is pos- 
sible to perform finite-size scaling in L to obtain proper- 
ties of the phase transition like critical exponents. We use 
a sparse-matrix techniqu o 32 ' 33 to enable the calculations 
of sizes up to L = 26. Note that using TM calculations 
one can in principle obtain the exact non-mean field ex- 
ponents of the phase transition by using sufficiently large 
L. This is unlike the CMF approach that we discussed 
in Sec. Mil where the exponents are mean field exponents 
(although for sufficiently large clusters the true exponents 
should emerge close to the transition, but such large clus- 
ters cannot be reached in practice). 

The critical exponents of the transition are obtained by 
calculating three types of scaled gaps based on the eigen- 
values of the TM. Each of them converges to a separate 
scaling dimension when system size L tends to infinity 
at the critical point. The antiferromagnetic scaled gap is 
defined as 

X h (T,g,L) = ^ln(^), (28) 

where Ao is the largest eigenvalue and Ai the largest 
eigenvalue in the subspace that breaks the symmetry of 
two neighboring sites, (i.e., corresponding to the stripe 
state), which means that the associated eigenvector V\ 
satisfies 

vi = -mil, (29) 

where R is the translation operator that translates the 
lattice by one unit along the axis of the cylinder. Thus, 
the scaled gap is also called the stripe scaled gap. For the 
system to host this stripe order, we have to restrict the 
system to even L. In addition, the eigenvector also bears 
odd parity when the system is reflected about the center, 
but is invariant under global spin flips. 



The two other scaled gaps are defined as 

X tl (T,g,L) = ^ln(^), (30) 

and 

X t2 (T,g,L) = ±ln(^j, (31) 

where A2 and A3 are the leading and sub-leading eigen- 
values associated to eigenvectors that are invariant under 
the lattice translation and global spin flips. However, it 
is not a priori clear which of these gaps corresponds to 
the thermal scaling dimension, and what the physical 
meaning is of the other gap. 

According to finite-size scaling theory 3 ^ and conformal 
invariance^ the gap Xi(T, g, L) in the vicinity of a crit- 
ical point scales as 

X,(T,g,L) = X l + a(T-T c )L yt + buL y " + ■ ■ ■ ,(32) 

where i indicates one of the three gaps (i = h, i = t\, or 
i = t-i), Vt is the leading thermal exponent, u the leading 
irrelevant field, and y u is the associated irrelevant expo- 
nent. The constants a, b are unknown (not universal). 

We calculate the scaled gap Xh(T, g, L) and then nu- 
merically solve for T C (L) using the following scaling equa- 
tion: 

X h (T,g,L)=X h (T,g,L-2). (33) 

The solution T C {L) converges to the critical point T c as 
L — > 00 in the following way: 

T C (L) = T C + a'uL^- yt + • • • , (34) 

where a' is an unknown constant. We thus determined 
the critical points of the 3\-Ji model for various g values 
to good accuracy, and the results are listed in Table Q] 
These T c values extracted from the TM approach agree 
very well with our MC results. 

The scaled gaps Xh, X tl and X t2 at the solutions T C {L) 
are calculated for a sequence of systems up to L = 26. 
Generally speaking, these gaps should converge to the 
corresponding scaling dimensions, respectively, in the fol- 
lowing way: 

Xi{L) = X l + b'L y « +•■■ , (35) 

where b' is an unknown constant. However, at the 4-state 
Potts point, the irrelevant exponent y u is zero, i.e., the 
corresponding field is marginally irrelevant, which leads 
to the following multiplicative logarithmic correction to 
scaling^ 

, N b\ &' 2 ln(lnL) , 

X i (L)=X i + ^ z + ^_l + ..., (36) 

with b'i, b' 2 unknown constants. 

From the scaling analysis of our MC data (Sec. IIVI) . 
we already know that g* « 0.67. Fitting Xi(L) according 
to Eq. (|3"5"|) for g > 0.67, we obtain the scaling dimensions 
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FIG. 13. (Color online) The scaling dimensions as functions 
of the Coulomb gas coupling qr. The curves for the AT model 
are theoretical predictions, while points are numerical results 
for the J1-J2 model with numerical TM-estimated gn. 



Xh, X tl , X t2 . The convergences of our data is not very 
good. This is because the irrelevant exponent y u has a 
small absolute value even away from g*. This is also the 
case for the AT model when K is close to 1 (the 4-state 
Potts critical point). For g = 0.67, the scaling dimensions 
are estimated by fitting Xi(L) to Eq. (f56"| . The results 
of such fits are listed in Table HI 

It is remarkable that, for all g in the region [0.67,5], 
the ratio of X t2 and X tl is always close to 4. Comparing 
with the CG formula Eq. (j4]) describing the AT univer- 
sality class, we thus identify X t2 as the thermal scaling 
dimension, and X tl a scaling dimension corresponding to 
the polarization scaling dimension X p of the AT model. 
Meanwhile the striped scaling dimensions are close to 1/8 
for all g, which corresponds to the magnetic scaling di- 
mension of the AT model. We further obtain gn for each 
g using Eq. ([4]) which expresses X p as a function of the 
CG coupling gR. These results are also listed in Table Q] 
We plot our numerical results X tl ,X t2 versus in Fig. 
[mi together with the CG predictions of X t and X p for 
the AT model. Thus, the TM calculations give a picture 
consistent with the MC simulations when we consider the 
ratio of the scaling dimensions X t2 /X tl . The TM calcu- 
lations, however, converge very slowly with L near the 
4-state Potts point, which makes the method unsuitable 
for the extraction of g* itself. 



VI. CONCLUSIONS 

We have shown that the thermal transitions from the 
striped ordered phase in the J1-J2 model for the range 
g £ [g*, 00) can be fully mapped to the continuous phase 
transitions of the well-known AT model. The special 
point g* « 0.67 corresponds to the 4-state Potts uni- 
versality class and for g — > 00 the transition approaches 
the standard Ising universality class. We have provided 



TABLE I. Best estimates for the critical properties obtained 
using the TM method for the J\-J% model at various g values. 
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1 1 IT 


Y, — n/<) 
-A-h — 7/ z 


y, — 9 _ 1 /,/ 
si-t 2 — ^ r/ v 




Qr 


0.67 


0.8335(2) 


0.12(1) 


0.50(2) 


0.12(1) 


4 


0.70 


0.7758(1) 


0.12(1) 


0.57(1) 


0.14(1) 


3.5 


0.75 


0.69866(5) 


0.12(1) 


0.630(5) 


0.155(5) 


3.17 


1.00 


0.48029(5) 


0.123(5) 


0.803(5) 


0.199(5) 


2.5 


2.00 


0.22468(3) 


0.125(2) 


0.953(2) 


0.238(1) 


2.1 


5.00 


0.088406(5) 


0.125(2) 


0.993(1) 


0.248(1) 


2.01 



a numerical mapping between the critical lines of the two 
models, based on matching universal properties; critical 
exponents as well as order-parameter histograms. 

Interestingly, the 4-state Potts model and the neigh- 
boring transitions in the AT model show a pseudo-first- 
order behavior on finite lattices, though these transitions 
are rigorously known to be continuous. The energy and 
order-parameter histograms show double-peak structures 
near T c , with the distance between the peaks decreas- 
ing slowly to zero as the system size is increased. This 
feature of the Potts point and its neighborhood conse- 
quently leads to similar effects in the J1-J2 model as well, 
in the vicinity of g*. This feature was misinterpreted 
as indicative of first-order transitions in some previous 
studies. The frustrated Ising model exhibits this pseudo- 
first-order behavior for g* < g < 1. 

We further showed that the MC data of the J1-J2 
model at g* can be scale-collapsed by using the critical 
exponents of the 4-state Potts model, if logarithmic scal- 
ing corrections known to exist at this point are properly 
taken into account. CMF and TM calculations were also 
used to understand aspects of the phase transition. CMF 
on 2 x 2 and 4x4 clusters predict a multicritical point at 
g* fts 0.66; very close to the exact location based on the 
MC calculations. However, since these calculations are 
mean-field in nature, the universality class of the contin- 
uous exponents cannot be determined directly with this 
approach. The TM calculations, carried out on cylinders 
of infinite length and finite width, lead to accurate (well 
converged) results for the transition temperature T c and 
also present an alternative method (to MC simulations) 
to calculate critical exponents when g is not too close 
to g*. However, it is difficult to reliably compute the 
location of g* with the TM method based on accessible 
cylinder widths, due to effects of the logarithmic correc- 
tions discussed above. 

An open issue requiring further investigation is to un- 
derstand why the transitions in the whole region (1/2, g*) 
(especially near g = 1/2, where an unusual first-order 
transition occurs at T = 0) are so weakly first-order (un- 
less they are of some more exotic continuous kind, which 
cannot be completely ruled out). 

The CMF calculations indicate that there may be a 
narrow region of first-order transitions for g < 1/2 (for 
the ferromagnetic-paramagnetic transition). It is not 
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clear whether this is an artifact of the small cluster size 
(with the first-order behavior obtaining for a 4 x 4 cluster 
but not for 2x2). Further large-scale MC simulations 
and TM calculations in this region are called for. 
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